Example of a first-order Neel to Valence-Bond-Solid transition in two-dimensions 
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We consider the S — 1/2 Heisenberg model with nearest-neighbor interaction J and an additional 
multi-spin interaction Q3 on the square lattice. The Q3 term consists of three bond-singlet projectors 
and is chosen to favor the formation of a valence-bond solid (VBS) where the valence bonds (singlet 
pairs) form a staggered pattern. The model exhibits a quantum phase transition from the Neel state 
to the VBS as a function of Q3 / J. We study the model using quantum Monte Carlo (stochastic series 
expansion) simulations. The Neel-VBS transition in this case is strongly first-order in nature, in 
contrast to similar previously studied models with continuous transitions into columnar VBS states. 
The qualitatively different transitions illustrate the important role of an emerging U(l) symmetry in 
the latter case, which is not possible in the present model due to the staggered VBS pattern (which 
does not allow local fluctuations necessary to rotate the local VBS order parameter). 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Mg, 75.40.Cx 



I. INTRODUCTION 

Even though quantum phase transitions (QPTs) occur 
strictly at zero temperature (T), they can strongly influ- 
ence the behaviour of a system at T > in their vicinity^ 
Exotic QPTs that do not have any analogs in classical 
critical phenomena may thus leave widely observable sig- 
natures. A very interesting example of such a QPT was 
proposed by Senthil et al.jz who reconsidered the transi- 
tion between Neel and valence-bond-solid (VBS) phases 
in two-dimensional (2D) square-lattice S = 1/2 anti- 
ferromagnetic spin systems i2j4 They argued that subtle 
Berry phase effects can lead to an unconventional "decon- 
fined" quantum-critical (DQC) transition that is gener- 
ically continuous, while a standard Landau- Ginzburg- 
Wilson analysis would instead predict a first-order transi- 
tion (except at fine-tuned multi-critical points), since the 
two phases break unrelated symmetries. At the DQC 
point, new fractionalized degrees of freedom, S = 1/2 
spinous, become liberated^ while they are confined in 
the conventional phases on either side of the transition. 

An important question raised by the field-theoretical 
work on DQC points^ is whether there are microscopic 
models (hamiltonians) where the existence of this uncon- 
ventional quantum-criticality can be explicitly demon- 
strated. Such models, apart from allowing unbiased tests 
of the validity of the concept, can lead to further theoret- 
ical developments. This is also important in view of the 
work reported by Kuklov et al. in Refs [H and [f| where 
an action claimed to realize the low-energy physics of the 
non-compact (NC) CP 1 field-theory, argued by Senthil 
et al. to describe the DQC point^ was studied using 
Monte Carlo simulations of 3D classical current models 
with easy-plane and SU (2) symmetry. The transitions 
in both cases were argued to be generically weakly first- 
order, in contradiction with the proposal of Ref.0. Other, 
similar studies have given different results, however 2& A 
full understanding of the NCCP 1 theory and its ability 
to capture the physics of the Neel-VBS transition in 2D 
antifcrromagnets is still lacking. 



One approach to the DQC problem is to identify and 
study microscopic quantum spin hamiltonians exhibit- 
ing the Neel-VBS transition. In a recently proposed 
class of S = 1/2 "J-Q" models a geometrically un- 
frustrated multi-spin interaction of strength Q (consist- 
ing of two or more bond-singlet projectors) was intro- 
duced which competes with the Heisenberg interaction 
J. Quantum Monte Carlo (QMC) calculations in this 
case are free from the sign problems^ hampering stud- 
ies of frustrated systems (which also exhibit Neel-VBS 
transitions^). Calculations in the ground state*^ and 
at finite temperature,— have demonstrated behaviors 
consistent with a continuous Neel-VBS transition of the 
DQC type. In particular, the dynamic exponent z = 1 
and the anomalous dimension associated with the spin 
correlations is relatively large; 77 m 0.35. There does, 
however, appear to be significant corrections to scaling at 
the critical point, possibly of the logarithmic typeJ^ Such 
corrections, which have also been observed by studying 
the response of the spin texture to impurities^ raise 
questions on the presence of marginally irrelevant opera- 
tors (or, possibly, conventional operators with scaling di- 
mension very close to 0) in corresponding field-theories. 
It has also been claimed that observed scaling correc- 
tions indicate an eventual first-order transition^ which, 
however, is unlikely, considering the fact that the scaling 
corrections are well behaved over a large range of system 
sizesJ^ Some scaling corrections are always expected, 
and it can be noted that logarithmic corrections have 
been found in related gauge theories with fermionsJ^ 

In light of the controversy regarding the nature of the 
Neel-VBS transition, and in order to gain further insights 
into the unusual physics of DQC points, it is useful to 
also study similar transitions that are definitely of first 
order. This is the topic of the present paper. Ideally, one 
would like to have a model in which the transition could 
be tuned from strongly to weakly first-order, through a 
tricritical point into a continuous transition. No such 
tunable model is known so far, however. In this work, we 
provide an example of a Neel-VBS transition on the 2D 
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FIG. 1. (Color online) a) The interaction term Q3 involving 
three bond-singlet projection operators (shown with thicker, 
red lines) on the square lattice. All terms related by lattice 
translations and rotations of the shown instance of a product 
of singlet projectors are included in the Hamiltonian. Exam- 
ples of VBS patterns: b) staggered, c) columnar, and d) pla- 
quette. The singlets preferentially form on the thicker (red) 
bonds. The arrows in c) indicate how a local resonance of 
a pair of bonds between horizontal and vertical orientation 
corresponds to a plaquette in d). Such resonances corre- 
spond to fluctuations of the VBS angle, which is (j> — nn/2 
(n — 0, 1, 2, 3) for the columnar state and (f>nn/2 + ir /4 for the 
plaquette state. In the DQC scenario^ they lead to an emer- 
gent U(l) symmetry, i.e., a continuous circular-symmetric (f> 
as the transition into the Neel state is approached. The Q3 
term favors the staggered- type VBS b), where such angular 
fluctuations should always be small. 



square lattice that is strongly first order. We introduce a 
J-Q3 model in which the Q3 interaction consists of three 
singlet projectors arranged in a staggered fashion on the 
square lattice, as illustrated in Fig. [T]Ja). This interac- 
tion induces a staggered VBS pattern in which resonating 
valence bonds leading locally to a different (degenerate) 
VBS pattern are not favored. If such fluctuations are 
present they can effectively rotate the coarse-grained an- 
gle of the VBS order-parameter (as explained further in 
the caption of Fig. QJ, which has been explicitly observed 
in the J-Q models studied previously ^ 10 ' 1 7 In the DQC 
scenario, they are directly responsible for the emergent 
U(l) symmetry of the VBS order parameter in the neigh- 
borhood of the transition^ The absence of this feature 
in the model studied here brings it clearly outside the 
framework of DQC points, and a numerical confirmation 
of a different type of transition is then, indirectly, an ad- 
ditional piece of evidence in favor of a consistent DQC 
scenario in which emergent U(l) symmetry and spinon 
deconfinement should go hand-in-hand with a continu- 
ous transition. 

We here use the stochastic series expansion (SSE) 
QMC method with operator-loop updates^ to study the 
nature of the Neel- VBS transition in the staggered J-Q3 
model. We perform simulations at a fixed aspect ratio 
of inverse temperature ft J — L, as done previously for 
the standard J-Q 2 model in Refs. [H, [lj, and [O- We 



study the finite-size scaling properties of various physical 
quantities and contrast them with what is observed at 
the previously studied putative continuous DQCs. 

The rest of the paper is organized in following way: In 
Sec II, we define the model more precisely and present 
the results for the staggered magnetization, the corre- 
sponding Binder cumulant, the spin stiffness, and the 
VBS order parameter. We also consider the probability 
distribution of the VBS order parameter and use it to 
explicitly demonstrate phase coexistence. In Sec III, we 
determine the location of the critical point by using the 
crossing of the energies of the Neel and the VBS phases 
in the meta-stable region near the transition. We state 
our conclusions and discuss future prospects in Sec IV. 



II. MODEL AND ORDER PARAMETERS 



We consider the following Hamiltonian 



H 



CijCklCn 



(1) 



(ijklmn) 



where Si refers to a S = 1/2 spin at site i on the 2D 
square lattice and Cij denotes the singlet pair projection 
operator, 



4 Si • Sj , 



(2) 



between two nearest neighbors i and j. The Q3 term 
(where, in the notation of Ref. [l(J the subscript on Q 
refers to the number of singlet projectors in the product) 
is chosen in the particular manner illustrated in Fig [Ha), 
to favor the formation of the kind of staggered VBS illus- 
trated in FigQJb). Like the columnar and plaquette VBS, 
the broken symmetry of the staggered VBS is Z4. How- 
ever, this type of VBS is very different from its columnar 
or plaquette counterparts, since no local ring-exchange 
of singlets on closed loops [e.g., as illustrated in FigQTc) 
for a simple two-bond resonance] is possible in the ideal 
staggered VBS. This makes it highly unlikely for the ex- 
isting fluctuations of this kind of VBS to be associated 
with an emergent U(l) symmetry, which is a key char- 
acteristic of the DQC transition^ We will confirm this 
with simulation results below. 

We have also studied an interaction similar to the six- 
spin Q3 term but with only two singlet projectors, on 
two pairs of sites separated by one lattice spacing and 
shifted one step with respect to each other as in Fig.QJa). 
This interaction is not sufficient for destroying the Neel 
order, however, unlike the original J-Q model with the 
two singlet projectors inside 2x2 plaqucttes. In the latter 
case the resulting VBS in the extreme case of J = is 
also quite weak,— while adding one more singlet projector 
(with the sets of three projectors arranged in columns) 
gives a much more robust VBS order ^ 

To study the Neel- VBS phase transition in the stag- 
gered J-Q3 model ([1]), we measure quantities that are 
sensitive to the Neel order and the VBS order respec- 
tively At a continuous quantum phase transition, these 
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FIG. 2. (Color online) The squared staggered magnetization 
{(m z s ) 2 ) shown for different system sizes at inverse tempera- 
ture (3J = L. The vertical line at {Qz/J) c = 1.1933 is the 
estimated L — > oo transition point from crossings of meta- 
stable energy branches (Fig[7J). 
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FIG. 3. (Color online) The Binder cummulant of the stag- 
gered magnetization shown for different system sizes at in- 
verse temperature /3J — L. Note that the minimum of the 
Binder cummulant is negative for L > 8 and diverges to — oo 
as L — > oo based on these sizes. 



quantities should scale with the system size L according 
to non-trivial critical exponents, while at a first-order 
transition one would expect very different exponents re- 
lated to the dimensionality of the system as well par- 
ticular signatures of coexisting phases at the transition 
point. These signatures should apply when the linear 
dimension of the system L > £, where £ is the finite 
correlation length at the transition. 



A. Neel order 

The magnetically ordered Neel phase breaks the SU(2) 
rotational symmetry of the interaction hamiltonian H 
and can be characterized by measuring ((to^) 2 ), where 
m z s denotes the z-component of staggered magnetization 
of the system: 



m~ 



N ^ 



S z (r) cos(Q -r) 



(3) 



with Q = (-7T, 7r) the wave- vector corresponding to the 
Neel phase and N = L 2 . This quantity is diagonal in 
the S z basis used and can be easily measured in the SSE 
simulations. We measure the squared quantity ((to;:) 2 }, 
which, due to the spin-rotational symmetry of the hamil- 
tonian, is 1/3 of the full squared staggered magnetiza- 
tion (to 2 ). We show the data for different system sizes 
at f3J = L near the phase transition in Fig [2] As the 
system size is increased, we observe a jump developing 
in ((m z ) 2 ) that becomes more abrupt and rapidly ap- 
proaches the infinite- volume estimate of the critical point 
{Qzl ' J)c = 1.1933(1) for this model (indicated by the ver- 
tical line in Fig[5]and other figures). The value of (Q3/ J) c 
was obtained from the crossing of metastable energies of 
the Neel and VBS phases of larger systems; see Fig [7] and 



later discussion in Sec. Mil This kind of behaviour of the 
Neel order-parameter is already very suggestive of a first- 
order transition. Data for systems larger than L = 14 are 
not shown here because of the extremely long tunneling 
times between the co-existing (as we will show below in 
Sec. Ill B| ) Neel and VBS phases for such sizes in our simu- 
lations close to the phase transition, which makes it very 
difficult to obtain reliable expectation values. 

A quantity that is very useful for distinguishing be- 
tween first-order and continuous phase transitions is the 
Binder cummulant U2, defined for an 0(3) order param- 
eter asi£ 



U 2 



I' 



(«) 4 ) 
3<(TOf) 2 >; 



(4) 



With the factors used here, U2 — >• 1 in the Neel phase 
and U2 — > in the magnetically disordered phase (VBS 
in this case) when L — > 00. For a continuous phase tran- 
sition, the U2 curves for different system sizes intersect at 
the critical point (for sufficiently large L) and the value of 
U2 at the intercept normally lies in the interval (0, 
This property of the Binder cummulant is often used to 
accurately determine the location of the critical point 
for continuous phase transitions. However, for a first or- 
der transition, the Binder cummulant behaves in a com- 
pletely different manner, that was explained phenomeno- 
logically for classical transitions by Vollmayr et al^ For 
systems exceeding a certain length L m i n ~ £, the curves 
show a minimum which becomes more pronounced as the 
system size increases. The minimum value of U2 — > — 00 
as L — > 00 because of phase co-existence, and the posi- 
tion of the minimum approaches the transition point in 
the thermodynamic limit. This behaviour has been ob- 
served in previous studies of classical first order phase 
transitions [for example, see Refs I20I and [2l| . Indeed, the 
Binder cummulant for the J-Q3 model, graphed in Fig [31 
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behaves in a similar manner and strongly points to a first 
order phase transition. The negative minimum is present 
in the U2 curves for L = 6 and above and becomes deeper 
and sharper as L increases. Its location approaches the 
estimated (Qa/J) c - The fact that a minimum in \J% is 
not at all present for L = 4 and is barely negative for 
L = 6 allows us to estimate that the typical length-scale 
(the spin correlation length) at the first-order transition 
is approximately in the range £ = 4^6. We also measure 
the second moment spin correlation length £ a defined as 



6. 



S(ir, 7r) 



L 



2ttV S(tt + ^,tt) 



(5) 



where S(q) refers to the spin structure factor at the cor- 
responding wavevector q. We obtain £ a ~ 2 close to the 
critical point in the magnetically disordered VBS phase. 
However, when the spin correlation length is small, £ a 
can differ from the true correlation length £ based on the 
asymptotic decay of the spin correlations in real space 
(which is difficult to extract reliably). 

An important quantity characterizing the Neel phase is 
the spin stiffness p s ; the second derivative of the ground 
state energy E(<f>) (per spin) in the presence of a twisted 
boundary condition (with rotation of a boundary row or 
column by an angle </>); 



Ps 



d 2 E{(/)) 



(6) 



This quantity is obtained in the SSE simulations in 
the standard way using the fluctuation of the winding 
numbersj 22 i 23 Again, the finite-size behaviour should be 
very different for a continous z = 1 phase transition and 
for a first order transition. In the former case, p s ~ l/L 
for large enough system sizes at the critical coupling. 
Thus, p s L curves for different system sizes intersect at 
the critical coupling for a z — 1 continous critical point, 
given L is sufficiently large. However, at a first order 
transition, since there is phase coexistence, p s is finite 
(non-zero) at the transition point when L — > 00 and thus 
p s L — > 00 as the system size is increased and no unique 
crossing point is obtained. Fig [4] shows the behavior of 
both p s and p s L. Intersection between p s L curves for 
different L can indeed be seen close to the transition 
point, but the value at the intersection points (for, e.g., 
L and L + 2) increases with L. The p s curves become 
sharper and approach a step function with a disconti- 
nuity at (Q3/ 1 J) c as L increases. Interestingly, the spin 
stiffness p s actually increases on approaching the transi- 
tion from the Neel phase before sharply dropping off in 
the VBS phase. 



B. VBS order 

To provide further evidence for the first order nature 
of the transition, and to obtain a bench-mark for the 
finite-size scaling, it is useful to also consider observables 
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FIG. 4. (Color online) The spin stiffness p s (top panel) and 
p s L (bottom panel) for different system sizes. 



sensitive to the VBS order. We define the following order 
parameters for that purpose: 

1 



r 



D 1 



N 
(Dl 



(7) 

(8) 
(9) 



where again Q = (7r,7r). The squared order parameter 
D 2 is zero in the Neel phase and becomes non-zero in 
the VBS phase in the thermodynamic limit. In Fig we 
show the behaviour of D 2 near the transition. The jump 
in D 2 gets sharper with increasing system size, similar 
to the staggered magnetization in Fig [51 and its location 
approaches {Qs/J) c - Thus, the solid order parameter 
also shows a first order transition. Note that the order 
parameter we use here is not spin-rotationally invariant. 
One can also define D x and D y using rotationally invari- 
ant bond operators such as S(r) • S(r + x). There may be 
differences in the finite-size scaling properties^ but for 
our purposes here the rotational invariance is not neces- 
sary (and the z-component definition of D 2 is easier to 
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FIG. 5. (Color online) Size dependence of the VBS order 
parameter D 2 for different system sizes. 



measure with the SSE method). 

A strong signature of a first order transition is phase 
co-existence at the transition point, unlike in the case of a 
continuous transition. We already showed evidence of co- 
existence of a Neel and non-magnetic phase in the form of 
a negative Binder cumulant of the staggered magnetiza- 
tion in Fig. [2] The coexistence can be observed in a more 
direct way using the VBS order parameter. Consider 
the joint probability distribution function P(D X , D y ). In 
the Neel phase, this function is peaked at (0,0). In the 
VBS phase, P{D X , D y ) is peaked at (0, ±D) and (±D, 0), 
where D is finite, reflecting the Z4 degeneracy of the VBS 
state in a finite lattice. Note that (D x ,D y ) is here de- 
fined as a single point obtained on the basis of an equal- 
time simultaneous measurement of D x and D y , i.e., these 
operators are not averaged over the imaginary-time di- 
mension in the simulations. The full distribution can still 
of course be accumulated over several imaginary times. 

In FigH we show P(D x ,D y ) for L = 12 and Q 3 /J = 
1.22, 1.23, 1.24 for (3J = L. The coexistence of the Neel 
and the VBS phase is evident from the presence of peaks 
at both (0, 0) and (±D, 0), (0, ±D) at Q 3 /J = 1.23, while 
at Q 3 /J = 1.22 (Q 3 /J = 1.24), the Neel (VBS) phase 
dominates. Also note the absence of any U(l) ring like 
feature in the distribution shown in Fig[6l This should be 
contrasted with similar measures of the distribution func- 
tion in the J-Q models with columnar VBS states close 
to the critical point, where the enlarged U(l) symmetry 
is very evident i & 10 i 17 In the original J-Q model with two 
singlet projectors, the VBS does not seem to get pinned 
to the four Z4 symmetric angles even at J = for the 
system sizes accessible.— In the modified J-Q model with 
three singlet projectors forming columns,— the change 
in the shape of the VBS order parameter distribution 
from U(l) close to the transition to Z4 deep in the VBS 
phase can be clearly observed, however. The enlarged 
U(l) symmetry arises in DQC theory due to the (dan- 
gerously) irrelevant Z4 symmetry-breaking term at the 



critical points Away from the critical point, the symme- 
try is only approximate but large system sizes are needed 
to observe that (i.e., L has to exceed the length scale A 
which is larger than the standard correlation length and 
determined by the dangerously irrelevant operator). 

Within the DQC framework, the approximate U(l) 
symmetry near the critical point can be thought of in the 
following manner £ the dangerously irrelevant Z4 pertur- 
bation only produces a small energy difference between 
the columnar and plaqucttc VBS which vanishes at the 
critical point, and gives rise to a corresponding large 
length scale slightly away from it within which the magni- 
tude of the VBS order parameter is formed, but its angle 
(which can be defined exactly as we did above in terms 
of D x and D y ) is not pinned in any particular direction. 
This feature has also been observed in U(l) symmetric 
spin models where either a very weak first order tran- 
sition^ or transition with unusual finite-size scaling^ 
takes place. However, when the VBS is staggered, there 
is no competing solid that is energetically close because 
of the absence of local ring-exchange moves of the sin- 
glets (or dimers) in the ideal staggered solid. This does 
not allow the emergence of an approximate U(l) symme- 
try near the transition, for which local fluctuations of the 
VBS order parameter are necessary, and puts it outside 
the framework of DQC points even though the phases 
have the same broken symmetries. Note that the value 
of D 2 in the VBS phase after the discontinuous jump 
(Fig [5]) is close (~ 74%) to that of an ideal staggered 
solid (D 2 = 0.015625), which motivates us to classify it 
as a strongly first-order transition. 

We should point out here that the order parameter 
distribution P(m z s ) of the Neel order parameter does not 
show a clear peak structure at coexistence, because of the 
spin-rotational averaging when measuring just one of the 
three components of the staggered magnetization. The 
distribution is not sharply peaked in the Neel state. In 
principle one could measure deviations from the rotation- 
ally averaged distribution expected^ for a single phase, 
but this signal is not as clear as the one seen above for the 
VBS order parameter. In the Binder cumulant this prob- 
lem is avoided, because only even powers of m s are used 
and they can be trivially related to the corresponding 
powers of m z s . In principle one could also measure the x 
and y components of the staggered structure factor, and 
compute the distribution P(\m s \), but since these are off- 
diagonal operators in the basis used there are ambiguities 
in how to define m x s and m y s for a given configuration in 
which ml is also measured. 



III. DETERMINATION OF THE TRANSITION 
POINT 

The location of a phase transition in the thermody- 
namic limit may be determined by using finite size scal- 
ing arguments. For first-order transitions at finite tem- 
perature, the peak position of the specific heat C v , sus- 
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FIG. 6. (Color online) The probability density P(D x ,D y ) 
shown for L = 12 at Q 3 /J = 1.22, 1.23, 1.24 and /3J = 12. 
The maxima present both at (0,0) and (±D, 0), (0, ±D) at 
Qs/J = 1.23 show phase coexistence of the Neel and VBS 
orders. 



ccptibility \ and the minimum of the Binder cummulant 
of the order parameter approach^! the thermodynamic 
transition temperature T C (L — > oo) as L~ d , where d is 
the dimensionality of the system. However, it is known 
from numerical studies of classical phase transitions (e.g 
in q > 4 Potts models in 2D) that large system sizes 
may be needed to observe this simple scaling of T C {L). 
For example, in Ref. H^, both the q = 8 and q — 10 Potts 



model (both of which exhibit first-order transitions) were 
studied for L < 50, and strong deviations from the L~ d 
scaling was observed for the q = 8 case at these system 
sizes. Continentino and Ferreira^ have proposed that for 
a quantum phase transition, similar definitions of g c (L) 
will approach g c {L — > oo) as L~( d+Z \ where g is the (non- 
thermal) parameter driving the quantum phase transi- 
tion and z is a dynamic exponent (this reflects the equiv- 
alence of the d-dimcnsional quantum model at T = to 
a (d + z)-dimensional classical model at finite tempera- 
ture). In our problem, the presence of Neel order at the 
critical point and using the geometry with (3 J = L imply 
that finite-size scaling is governed by low-energy linearly 
dispersing spin waves and therefore z = 1. The loca- 
tion of the minimum of the Binder cummulant XJi can be 
taken as a possible definition for such a psuedo-critical 
coupling g c {L). With this definition, we have not been 
able to check the above mentioned scaling for our model 
explicitly because of limited system sizes, i.e., the data 
do not follow pure power laws. To obtain reliable data 
for bigger system sizes near the critical point, where the 
tunneling time between the competing phases becomes 
very large, will possibly require extended ensemble meth- 
ods (which have been developed for quantum system in 
Refs. l3(il - [32l ), which we leave for future work. 



Interestingly, the long tunneling times associated with 
strong first order transitions can actually be used to lo- 
cate the value of the critical point accurately. This can 
be done by determining the average energy of the two 
different phases^ in the following manner: To determine 
the energy of the Neel phase, we start the simulation at 
a value of Q3/ J inside that phase and then keep increas- 
ing Qzj J to go to the VBS phase, while measuring the 
average energy at each such Q3/J. In this sequence, we 
always use the last configuration generated in the run 
at a particular Q3/J as the starting seed for the next 
value of Q-ij J ■ Given the system is sufficiently big so 
that the tunneling time is large compared to the simula- 
tion time, this procedure ensures that the system stays in 
the metastable Neel phase close to the critical point even 
when it has been crossed. We then repeat this process in 
the reverse direction starting from a value of Q3/J inside 
the VBS phase and decreasing the value of Q3/J. The 
crossing point of the two energy branches then gives the 
location of a psuedo-critical coupling (Q 3 / J) C (L) at that 
system size. In Fig [7] (top panel), we show the results 
of this procedure for L = 32 and L = 36. The curves 
can be fitted very well to quadratic polynomials (the 
small quadratic part is added to improve x 2 / dof) and 
the crossing point can be accurately determined. These 
crossing points are indeed consistent with a scaling of the 
form (Q 3 /J) C (L) = (Q 3 / J) c — A/L 3 for the larger system 
sizes studied (L = 24 — 36) (see Fig [7] (bottom panel)) 
assuming that z = 1. We obtain (Qs/J),- = 1.1933(1) 
using this extrapolation. 
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FIG. 7. (Color online) Critical coupling (Qz/ J) c (Li) extracted 
from the crossing of the energy branches (top panel) of the 
Neel phase (labeled as fw) and the VBS phase (labeled as rev) 
for a system of size L = 36 (open symbols) at ft J = 36 and 
L = 32 at /3J — 32 (closed symbols). The bottom panel shows 
a fit to the form (Q 3 /J) C (L) = {Qz/J) c + A/L 3 to extract the 
thermodynamic value of {Qz/J)c = 1.1933(1). 



IV. CONCLUSION 



sociatcd ring-like features in the probability distribution 
of the VBS order parameter close to the critical point. 
We also use the crossing of the energy branches of the 
two phases for larger system sizes to estimate the ther- 
modynamic value of the critical point. We observe that 
these crossing points (for L = 24 — 36) indeed converge 
to the thermodynamic value of the critical point as 1/L 3 , 
expected from finite-size scaling for first order transitions 
with /3 r~j L. It will be interesting to test similar predic- 
tions of finite-size scaling theory of first-order quantum 
phase transitions^ which should depend solely on the di- 
mensionality of the system d and the dynamical exponent 
z, by measuring the order parameters and the Binder 
cummulant using quantum Monte Carlo simulations in 
combination with extended ensemble methods 3 ^— for 
bigger systems that may avoid the large tunneling times 
between the two co-existing phases near the transition. 
It will also be interesting to see whether there is a finite- 
temperature crossover to a continuous phase transition 
between a VBS and a disordered phase. Given the strong 
first order nature of the transition, it will likely sur- 
vive at low T as well. Though the transition is strongly 
first-order, there are still fluctuations present near the 
critical point, e.g. D 2 is 74% of its maximal value in 
the VBS phase in the vicinity of the transition and the 
length scale of appeareance of first order signatures in 
the Binder cummulant of the staggered magnetization is 
Lwi Mean-field theorie o 34 ' 35 have not been success- 
ful for the original J-Q model. It will be interesting to 
see whether a reasonable description of this transition 
and its location can be obtained using similar techniques 
(e.g. mean field theory based on bond-operator formu- 
lation for S = 1/2 introduced by Sachdev and Bhatt in 
Ref [36h and this is presently being investigated. Finally, 
none of the finite-size signatures of first-order transitions 
like negative divergence of the Binder cummulant, cross- 
ing of energy branches etc presented in this study, have 
been observed in previous simulations of SU(2) symmet- 
ric J-Q type models with a columnar VBS phase even for 
very large system sizes (L up to 256).— 



In conclusion, we have considered a modified "J-Q3" 
type Hamiltonian, where the Q3 term is chosen to sta- 
bilize a staggered VBS solid. We demonstrate the pres- 
ence of a strong first-order quantum critical point that 
seperates the Neel and the VBS phases. In the VBS 
phase, there is no approximate U(l) symmetry and as- 
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